Geometric Integration of ODEs Using Multiple Quadratic Auxiliary Variables

نویسندگان

چکیده

We present a novel numerical method for solving ODEs while preserving polynomial first integrals. The is based on introducing multiple quadratic auxiliary variables to reformulate the ODE as an equivalent but higher-dimensional with only integrals which midpoint rule applied. can subsequently be eliminated yielding midpoint-like original phase space. resulting shown discrete gradient method. Furthermore, averaged vector field obtained special case of proposed extended higher-order through composition and illustrated number examples.

برای دانلود باید عضویت طلایی داشته باشید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Geometric integration of conservative polynomial ODEs ✩

We consider systems whose Hamiltonian is of the form H(q,p) = 2p + V (q), where the potential V is either cubic or quartic with no cubic terms. For most of these systems (in the measure sense) we give an explicit numerical integration method that preserves both phase space volume and the value of the Hamiltonian. This is exemplified in the Hénon–Heiles system. An application is to the hybrid Mo...

متن کامل

Exterior algebra, ODEs and geometric integration References

References [1] A. Adrover, S. Cerbelli & M. Giona. Exterior algebra-based algorithms to estimate Liapunov spectra and stretching statistics in high-dimensional and distributed systems, Int. [4] L. Allen & T.J. Bridges. Flow past a swept wing with a compliant surface: stabilizing the attachment line boundary layer, Stud. Hydrodynamic stability of the Ekman boundary layer including interaction wi...

متن کامل

Survival analysis using auxiliary variables via non-parametric multiple imputation.

We develop an approach, based on multiple imputation, that estimates the marginal survival distribution in survival analysis using auxiliary variables to recover information for censored observations. To conduct the imputation, we use two working survival models to define a nearest neighbour imputing risk set. One model is for the event times and the other for the censoring times. Based on the ...

متن کامل

Membrane geometry with auxiliary variables and quadratic constraints

Consider a surface described by a Hamiltonian which depends only on the metric and extrinsic curvature induced on the surface. The metric and the curvature, along with the basis vectors which connect them to the embedding functions defining the surface, are introduced as auxiliary variables by adding appropriate constraints, all of them quadratic. The response of the Hamiltonian to a deformatio...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

ژورنال

عنوان ژورنال: SIAM Journal on Scientific Computing

سال: 2022

ISSN: ['1095-7197', '1064-8275']

DOI: https://doi.org/10.1137/21m1442644